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TOTAL AND JET BLANDFORD-ZNAJEK POWER IN PRESENCE OF ACCRETION DISK 

Jonathan C. McKinney 1 

July 14, 2005 

ABSTRACT 

A rotating black hole probably powers ultrarelativistic jets in gamma-ray bursts, relativistic jets from some 
active galactic nuclei, and jets from some black hole x-ray binaries. Prior estimates of the power output of 
a black hole have assumed an infinitely thin disk, a magnetic field based upon a slowly rotating black hole, 
and have not self-consistently determined the geometry or magnitude of the magnetic field for a realistic ac- 
cretion disk. We provide useful formulae for the total and jet Blandford-Znajek (BZ) power and efficiency 
as determined self-consistently from general relativistic magnetohydrodynamic numerical models. Of all jet 
mechanisms, we suggest that only the BZ mechanism is able to produce an ultrarelativistic jet. 

Subject headings: accretion disks, black hole physics, galaxies: jets, gamma rays: bursts, X-rays : bursts 
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1. INTRODUCTION 

For Poynting-dominated jets, where field lines tie the black 
hole to large distances, the energ y flux is determined by the 
Blandford-Znajek (BZ) proce ss (Blandford& Znaiek[^)7J|) 
(for a review see iRuffini & Wilson! 1 19751 iRees et a flll982t 
Begelman et alJ Il984 iBeskin & K uznetsova| j200(A iPunslv 
2001; Mc Kinn ev & Gammiel2004lLevinsorJ2005l) . 

The BZ effect depends on the magnetic field strength 
near the black hole and the Kerr black hole spin parame- 
ter a/M, where — 1 < a/M < 1. Self-consistent production 
of a relativistic Poynting jet likely requires a rotating black 
hole accreting a thick disk w ith a disk height (H) to ra- 
dius (R) ratio of H/R > 0.1 llGhosh & Abramowiczlll997t 
iLivio. Qgilvie. & Pringlelll999t iMeierJ 120011) . As discussed 
below, a rapidly rotating (a/M ~ 0.5-0.95) black hole ac- 
creting a thick (H /R > 0.1) disk is probably common for jet 
systems. 

However, prior estimates of the BZ power output assume 
the presence of an infinitely thin (H/R ~ 0) disk, only ap- 
ply for a/M ~ 0, and do not self-consistently determine the 
magnetic field strength or field geometry. The purpose of 
this paper is to provide useful formulae for the total and jet 
Blandford-Znajek power for arbitrarily rapidly rotating black 
holes accreting a realistic disk. We also discuss the domi- 
nance of the BZ effect over other relativistic jet mechanisms. 

2. BLACK HOLE ACCRETION SYSTEMS 



The accretion 
ing bl ack hole ' 



of a thick disk around a rapidly rotat- 
is often invoked as the engine to power 
GRBs (INaravan et alJI 1 992t IWooslevll 1 9931 i Paczynskill 1 99 
MacFadven & Wooslev 1999; Naravan et al. 2001; Broderic 
12004 . Typical GRB models invoke a rel atively thick (H/R ~ 
0.1- 9) disk jMacFadyen & Wooslevlll999t IPooham et alJ 
119991 iKohri et al.1 120051) . During the GRB event, the 
black hole forms with a/M ~ 0.5-0.75 an d evolves to a 
rapidly rotating state with a/M ~ 0.9 — 0.95 (INaravan et al.l 
1 19921 [M acFadven & Wooslev 1999t lShapiro & Shibatal2002t 
Shibata & Shapiro 2002). GRB models based upon internal 
shocks require an ultrarelativistic jet with a typical Lorentz 
factor of r ~ 100- 1000 in order to overcome the com- 
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pactness problem tLithwick & Sarill2001t |pkail 2005). while 
Compton drag models requ i re r ~ 20- 100 ( Ghiselli ni et alJ 
2000; Lazzati et al.1 12004 iBroderickl l2004lK Direct ob- 
servations of GRB afte rglow show evidence for relativis- 
tic motion (iGoodmanl 1 1 9971 iTavlor et alJ l2004albl) . Large 
Lorentz factors require a relatively large jet energy flux, 
which could be BZ-driven and Poynting-dominated rather 
than neutrino-a nnihilatio n-driven and enthalpy-dom inated 
(Meszaros & R eesll997tlP opham e t alJl999tlDi Matteo et ail 
2002; McKinnev 2005a b). Core-collapse of a rapidly rotat- 
ing star leads to an inner disk with a strong uniform (perhaps 
net) poloidal field. 

The accretion of a relatively thick (H/R ~ 0.9) disk around 
a rapidly rotating black hole is probably the engine that pow- 
ers jets from AGN and some black hole x-ray binaries. Both 
radio-loud AGN and x-ray binaries in the low-hard state show 
a correlation between radio and x-ray emission, which is con- 
sistent with radio synchrotron emission and hard x-ray emis- 
sion generated from Comptonization through a thick disk 
(Merl oni et aTJl2003l) . This suggests the disk is geometrically 
thick when a system produce s a jet, w here the disk is probably 
ADAF-like with H/R ~ 0.9 (iNaravan & Yill995l) . 

Based on Soltan-type arguments, AGN each probably 
harbor a rapidly rotating (a/M ~ 0.9-0.95) blac k hole 
JUrrv & Padovanil 119951 lElvis. Risaliti. & Zamoranil 120021 
iGammie. Shapiro. & McKinnevll2004 IShapiroll2005l). AGN 
are o b served to have jet s with V < 10 (lUrry $ Padoyanil 
1 1995 1 : iBiretta et ail 119991) even T ~ 30 (iBegelman et alJ 
1 1994: iGhisellini & Celottill2001t Uorstad et all 120011) while 
some observations impl y T < 200 (IG hisellini et alJ 119931 
iKrawczvnski et alJl2002tlKonopelko et alJ l2003). For exam- 
ple, the jet in M87 shows a large-scale opening angle of 10° 
with r ~ 6 (lJunor et alll999tlBiretta et alJ2002l) . AGN prob- 
ably a ccrete a uniform field from solar wind captu re or the 
ISM ( INaravan et all2003tlUzdenskv & Spruii2005l) . 

Bla ck hole x-ray binaries might hav e a/M ~ 0.5-0.95 
(Gammie. Shapiro. & McKinnev! 12 0041). while some may 
have a/M < 0.5 (Gierlihski & Done 200 #. X-ray binary sys- 
tems produce outflows an d jets (Mirabel & RodrigueJ ll999t 
McClintock & Remillard 2003). For example, black hole 
x-ray binary GRS 1915+105 has a je t with apparently su- 
perlum inal motion with V ~ 1.5 — 3 (M irabel & Rodriguez 
1 19941 iMirabel & Rodriguez! 119991: iFender & Bellonil 12004 
Kaiser et alT l2004l) . Solar- wind capture x-ray bi naries prob- 
ably accrete a uniform field ( Narav an" et all20 03). 



2 




0.1 r 

10" 2 t- 

= A A a A 
P 10- 3 t- A 



10" 4 
10- 5 
0.2 

0.15 



0.1 



0.05 



10 2 

% 10 

-5 
o 

1 

0.1 



total 



polar jets 



H — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1- 
E^uiparUtion Field 



0.2 



0.4 O.i 

a/M 



0.E 



FIG. 1. — Top panel: total (triangles: data, dashed line: fit) and jet (squares: 
data, solid line: fit) efficiency. Open points represent negative efficiencies. 
Middle panel: coefficient (r)) least squares fit to rj formulae. Bottom panel: 
normalized field (in Gauss) squared. 



3. THE BLANDFORD-ZNAJEK EFFECT 

Most authors est i mate the BZ power based upon the 
IB landlord & Znaiekl i ll 9771) model of a slowly spinning black 
hole threaded by a monopole-based magnetic field and accret- 
ing an infinitely thin disk, which gives 



Pbzm ~ Po(B r [G]T (n 2 H / c )A, 



(1) 



where B r is the radial field strength, r g = GM/c 2 , f2# = 
ac/ilMrn) is the rotation frequency of the hole, r# = r g (l + 
y/l — (a/M) 2 ) is the radius of the horizon for angular mo- 
mentum J = aGM/c, and the dimensionless Kerr parame- 
ter has -1 < a/M < 1. The parameter P = 0.01 - 0.1, 
where the uncertainty in Pq arises because the strength of the 
magn etic field is not self-con s istently d etermine d (see, e.g., 
MacDonal d~& Thornd Il982bt iThorne & MacDonaH 119821: 
Thorne, Price, & MacDonald 1986JK F orce-free numerica l 
models agree with the above BZ model (Komissarov I200I . 
GRMHD numerical models of slowly spinning, accret- 
ing black holes mostly agree with the BZ model for the 
nearly force-free funnel regio n of the Poynting-dominated jet 
(McKi nnev & Gammiel20 04). 

The force-free solution for the monopole BZ flux is oc 
sin 2 9 iBlandford & Z naiek 1977), but the accretion of a rel- 
atively thi ck disk di minishes total BZ power output sub- 
stantially (McKinnev & Gammie 2004). This is because 
the electromagnetic energy accreted as a disk dominates the 
energy extracted. Some black hole spin energy does es- 
cape into the diffuse part of the corona, so the coronal 
outflo w has more Poynting flux for faster sp inning holes 
JMcKinnev & Gammiel20 04: Kroli k et al.l2005l) . For rapidly 
rotating black holes, the field is no longer monopolar and a 
significant amount of flux is generated closer to the nearly 
force-free poles. 

HARM ( Gammie et al. 2003a|) was used to evolve a se- 
ries of otherwise identical GRMHD models with spin a/M = 



-0.999-0.999 and H/R = 0.1,0.2,0.5,0.9 to determine 
the total BZ power (P tot ), jet BZ power (Pj e t), and field 
strength (B r ). A similar ser ies of models were studied in 
iMcKinnev & Gammid (12004), and they give a description of 
the model setup, limitations, and related results. The evolved 
field geometry is relevant to most black hole systems and cor- 
responds to a turbulent disk field with a self-consistently gen- 
erated large-scale flux threading the black hole. 

Figure[0shows the data and fits described below for H/R = 
0.2. There is a weak dependence on H/R> 0.1 for the jet BZ 
power or efficiency since the thicker the disk, the less solid 
angle available to the jet, but the field strength there is larger 
in compensation. The non-jet results for other H/R will be 
presented in a separate paper. 

The total (disk+corona+jet) BZ power efficiency in terms 
of the mass accretion rate (M) for a/M > 0.5 is well fit by 



Vtot = 



Mc 2 



14.8 



Q H [a/M=l] 



(2) 



where the coefficient (ij) is obtained by a least squares fit. 
Net electromagnetic energy is accreted for a/M < 0.4 (in- 
cluding retrograde) when an accretion disk is presen t. Thi s 
fact surprisingly agrees with the thin disk study of Li (2000), 
who find that r\ tot > only if a/M > 0.36, corresponding to 
fin > fijeflSCO], the Keplerian angular velocity at the inner- 
most stable circula r orbit (ISCO). The sparse spin study of 
iKrolik et all (120051) is in basic agreement with our r\ tot . No- 
tice that the spin dependence of our r\ to t is consistent with 
IMcKinnev & Gammid J2004I) for their fit given by equation 
61 of data shown in figure 11, where they found r\ toi oc (2 — 
r H ) 2 oc (a/M) 4 . Their fit coefficient of 6.8% was inaccurate 
for a/M fts 1 since they had no model beyond a/M = 0.97 and 
they included points with a/M < 0.5 that do not fit well. 

The nearly force-free jet region contains field lines that tie 
the black hole (not the disk) to large distances, leading to the 
BZ-effect and a jet. For a/M > 0.5, 



(3) 

n H [a/M=\} / ' 

over both polar jets. If a/M w 0.9, then w 1% of the accreted 
rest-mass energy is emitted back as a Poynting jet. 

The horizon value of B r = F , where F is dual 
of the Faraday, determines the black hole power output 
jMcKinnev & Gammiel2004l) . For all a/M > 0, the total and 
jet fields are 

4 

(4) 



™^ . 6 + 20 

POJiskC 1 

(B r jet [G]) 2 

PO.diskC 2 



1.0 + 81 




(5) 



where the equipartition field satisfies (B''[G]) 2 /(87r) 
Po.diskC 2 , where p Q j isk = Mt g / r 3 g and t g = GM/c 3 . Hence, 

P tot «7.4 x 10- 3 ((B, r o ,[G]) 2 r 2 c-0.6Mc 2 ), 



Pjet'' 



;8.4 x 10- 4 ((fi; Vf [G]) 2 r 2 c- l.OMc 2 ), 



(6) 



(7) 



where since the field is determined self-consistently, no ex- 
plicit spin dependence appears. This demonstrates the com- 
petition between electromagnetic energy extraction and accre- 
tion. Notice that no direct comparison can be cleanly made to 
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FIG. 2. — Field type #9 dominates in GRMHD numerical models and is 
associated with Blandford-Znajek effect. Types #1,2,3,5,6 dynamically im- 
portant. Type #4 transient. Types #7,8 not dynamically stable. 



equation[2due to the presence of two ambiguities Po and B r , 
while our formulae have no ambiguities. We suggest using 
the power output as given by equations [2] and [5] The black 
hole mass accretion rate must be determined independent of 
a /M using a model-dependent study. 

The coefficient in the formulae above depends on the 
type of accreted field geometry. As an extreme exam- 
ple, the accretion of a net vertical field leads to an in- 
creas e in the net electromagnetic efficiency by a factor of 
five jMcKinnev & Gammid 12004). Also, the accretion of 
a net toroidal field lea ds to negligible energy extraction 
l lDe Villiers et al . 2005a). Future studies should focus on the 
physical relevance, stability, and long-temporal evolution of 
accreting net toroidal and vertical fields with realistic pertur- 
bations rather than exact symmetries. Accretion of a net ver- 
tical field has been studied with nonrelativistic MHD simula- 
tions ( Igumens hchev et alJ2003l) . 

4. BZ JET POWER FOR COLLAPSAR MODEL 

For the collapsar model with black hole mass M ~ 3M 
feeding at an accretion rate of M = OAMq/s, the magnetic 
field at the poles of an a/M ~ 0.9 black hole is B r ]et w 10 16 G 
and po,disk ~ 3.4 x 10 10 gcm" 3 . This gives a per polar axis 
Poynting flux jet energy of Pj et rj 10 51 erg s" 1 . Notice that the 
neutrino annihilation jet luminosity f or such a collapsar mo del 
gives Luu^njet ~ 10 5() - 10 51 erg s" 1 (iPopham et all 19991) . so 
these processes are likely both important. However, real- 
istic models suggest that Poynting flux domi nates neutrino- 
annihilation energy flux (McKinnev 2005a b). Similar esti- 
mates can be made for AGN and x-ray binaries. 

5. DOMINANCE OF BLANDFORD-ZNAJEK EFFECT 

Figure |2] shows the possible types of field geomet ries in 
the disk (see also, e.g.. lBlandfordll2002l Iffirose et al.l l2004). 
Field type #1 correspon ds to the Balbus-Hawley instability 
( Barb us&Hawlevl 19911) . which is present in our simulations. 

Field type #2 corresponds to models for which the field 
ties material inside the ISCO to the outer disk (Gammie 
fl999t lKrohki rr999). As they predicted, unlike in the a- 
viscosity model, there is no f eature at the ISCO or a di- 
rect plunge into th e black hole (M cKinnev & Gammiell2004l 
iKrolik et alJ 120051) . which impacts any radiative model of 
the inner-radial accretion disk. Field type #3 corresponds 



to models that conside r the role of th e black hole a nd disk 
on the disk efficiency JGammielll999l lKroliklll999l) . They 
suggested that efficiencies of order unity or higher could be 
achieved by extracting energy from the hole. This field type 
is pres ent, but the disk efficiency i s near the thin disk effi- 
ciency (McKinnev & Gammie 2004). Thus, surprisingly, the 
magnetic field and disk thickness play little role in modifying 
the disk efficiency. In contrast, the angular momentum ac- 
creted is re duced by magnetic f ield lines that tie from the dis k 
to the hole dMcKinnev & Gammiel2004tlKrolik et alJ2005l) . 

Field types #4 and #5 correspond to surface reconnections. 
Type #4 geometries are temporary and type #5 are common. 
Thus, reconnection efficiently removes large loops that tie the 
disk to itself. Field type #6 corresponds to the Blandford- 
Payne type model (Blandford & Pavne 1982). We find that 
the lab-frame \B\ oc r" 5 / 4 as in their model, but the lab-frame 
p ex r°'° instead of their p oc r" 3 / 2 and there are few such field 
lines present in a stable configuration due to the inner-radial 
corona being convectively unstable and magnetically unsta- 
ble to magnetic buoyancy. Thin disks likely have more stable 
surfaces that might allow for a stable wind. 

Field type #7 correspond s to coronal out flows or 
ergospheric-driven winds (Punslv & Coronitil 1990 a b). There 
are no dynamically stable field lines that tie the inner-radial 
disk to large distances. Even for a/M = 0.999, no additional 
Poynting flux is created in the ergosphere and the electromag- 
netic energy at infinity completely dominates the hydrody- 
namic energy at infinity associated with the MHD Penrose 
proces s, in basic agreement with the r esults of iKomissarovl 
(120051) and c ou nter to the r esults of iKoide et all (120021) : 
IPunslvl (12005b lKoide etal.1 ( l2002h evolved for much too 
short a time. IPunslvl (120051) used the 3D GRMHD near- 
horizon results of lKroliketa but their near-horizon 
results could have numerical artifacts associated with their 
use of Boyer-Lindquist coordinates. However, there is a con- 
vectively and magnetically unstable, self-consistently mass- 
loaded, colli mated, mildly relativistic (v / c < 0-95) coro- 
nal ou tflow dMcKinnev & Gammiel 120041: iDe Villiers et alJ 
12005 al) . A rotating black hole is not requ ired for nonrelativis- 
tic (v /c < 0.6) coronal outflows ( M cKinnev & Gammiel2002l 
2004). 

Field type #8 corresponds to Uzdenskv (2005) type mod- 
els. These field geometries appear rarely and do not transfer 
a significant amount of energy or angular momentum. Such 
geometries may be more important for thin disks. 

Field type #9, the dominant feature, is associated with the 
Blandford-Znajek model. Since the magnetic field confines 
the disk matter away from the polar region, the rest-mass flux 
there is arbitrarily low. The large BZ flux to low rest-mass 
flux ratio can translate into an arbitrari ly fast jet. The mass- 
loading of this jet is considered in McKinnev ( 2005a Q). 

Notice that for accretion models with a net verti- 
cal field, the resulting stru cture is essentially identical 
( McKinnev & Gammie 2004). Reconnection efficiently 
erases the initial geometrical differences. 

6. CONCLUSIONS 

Typical BZ power output estimates assume an infinitely thin 
disk, a slowly rotating black hole, and do not self-consistently 
determine the magnitude or geometry of the magnetic field. 
We use GRMHD numerical models to self-consistently deter- 
mine the total and jet BZ efficiency when a disk is present. 
There is a significantly stronger dependence on black hole 
spin than prior estimates suggest. 
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Near the rotating black hole, the field geometry of the ac- 
cretion system is dominated by the field that leads to the BZ- 
effect. Since the polar region is magnetically confined against 
disk material and the BZ power is large, the jet Lorentz factor 
can be arbitrarily large. This is unlike disk-related jet mecha- 
nisms that are directly loaded by disk material. 
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